Probing the nature of dark matter with Cosmic X-rays: 
Constraints from "Dark blobs" and grating spectra of galaxy clusters 
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Gravitational lensing observations of galaxy clusters have identified dark matter "blobs" with 
remarkably low baryonic content. We use such a system to probe the particle nature of dark matter 
with X-ray observations. We also study high resolution X-ray grating spectra of a cluster of galaxies. 
From these grating spectra we improve the conservative constraints on a particular dark matter 
candidate, the sterile neutrino, by more than one order of magnitude. Based on these conservative 
constraints obtained from Cosmic X-ray observations alone, the low mass (m 3 < 10 keV) and low 
mixing angle (sin 2 (28) < 10~ 6 ) sterile neutrino is still a viable dark matter candidate. 

PACS numbers: 



I. INTRODUCTION 

The cosmological dark matter abundance is firmly es- 
tablished through observations of the cosmic microwave 
background and of the large scale structure of the uni- 
verse [H . This is complemented by measurements of 
dark matter on smaller scales by studies of, e.g. the ro- 
tation curves of galaxies, gravitational lensing by galax- 
ies and clusters of galaxies , the velocity dispersion of 
galaxies in clusters of galaxies, and X-ray emitting hot 
gas in clusters of galaxies However, the particle na- 
ture of dark matter remains a puzzle. 

There are numerous dark matter candidates, among 
which the sterile neutrino is particularly well motivated. 
The sterile neutrino is a natural dark matter candidate in 
a minimally extended standard model of particle physics 
[Tl| and it provides solutions to other problems: the 
masses of the active neutrinos 0, the baryon asymme- 
try of the Univer se |^| , and the observed peculiar veloci- 
ties of pulsars |3, Il0( . Sterile neutrinos participate in the 
flavour-mass eigenstate oscillations of the active standard 
model neutrinos, and are thereby allowed to decay ra- 
diatively through a two-body decay with photon energy 
predicted to lie in the X-ray range (E^ — m s /2, where 
m s is the rest mass of the sterile neutrino) . This renders 
it a testable dark matter candidate [l2|. 

The decay rate of any dark matter candidate with a 
radiative two-body decay can be constrained from ob- 
servations of dark matter concentrations (for references 
see The strongest constraints are obtained 

from studying dark matter dominated regions, and with 
instruments with high spectral resolution (since the de- 
cay line is expected only to suffer negligible broadening) . 

Recent gravitational lensing observations of the mass 
distribution in galaxy clusters p| have identified cluster 
scale dark matter "blobs" with very low baryonic content. 
This allows for the novel possibility of using such almost 
pure dark matter blobs to probe the particle nature of the 
dark matter Below we analyze X-ray observations of 



the dark matter blob in the cluster of galaxies Abell 520. 
Also, we analyze high resolution X-ray grating spectra of 
the cluster of galaxies Abell 1835, allowing us to improve 
constraints from earlier studies by more than one order 
of magnitude. 



II. X-RAY DATA ANALYSIS 



When a spectrum has been obtained from an observa- 
tion of a given dark matter dense region, there are differ- 
ent ways of searching for a hypothetical mono-energetic 
emission line and to determine an upper limit on the flux 
from decaying dark matter particles. The simplest and 
most conservative method is the "slice method," where 
the energy range of the spectrum is divided into bins of a 
width equal to the instrumental energy resolution (2a), 
and all of the X-ray flux in a particular bin is determined 
from a model fitted to the spectrum. The slice method 
is very robust, as the physics behind the fitted model 
is irrelevant, and the method does not require any as- 
sumptions about the X-ray background, but regards all 
received flux as an upper limit for the flux from decay- 
ing dark matter. This is despite the fact that the total 
flux is known to consist of several contributions; the cos- 
mic X-ray background from unresolved sources, the X- 
ray emission from the intra cluster medium, the Milky 
Way halo, and the instrumental background. The "slice 
method" takes into account that an emission line from 
decaying dark matter could "hide" under a line feature 
in the spectrum Il6|. Other methods for constraining the 
flux (for example [T3) can give stronger, but less robust, 
results (for a discussion of different methods see 0]). In 
this study, we conservatively use the slice method. 



2 




III. GRATING OBSERVATIONS OF THE 
GALAXY CLUSTER ABELL 1835 



A good spectral resolution is required to search for a 
mono-energetic emission line. The high spatial resolu- 
tion of the Chandra X-ray telescope can be turned into a 
very high spectral resolution (» 5eV) by deflecting the 
incoming photons in a grating, as the deflection angle is 
highly sensitive to the photon energy. The ACIS-HETG 
instrument consists of two gratings, which on demand 
can be placed between the mirrors and the ACIS CCDs: 
the High Energy Grating (0.8-9.0 keV), HEG, and the 
Medium Energy Grating, MEG (0.4-5.0 keV) 33]. 

When the incoming photons are deflected in a grating, 
the information of their spatial origin is lost. This makes 
it impossible to optimize the ratio of expected dark mat- 
ter signal to noise from X-ray emitting baryons in the 
observational field of view (as described by [la, 13)- ^ n 
this study, the cluster of galaxies Abell 1835 was targeted 
because most of the emission from the cluster comes from 
a region close to the observational axis making grating 
spectroscopy possible. Abell 1835 has a luminosity dis- 
tance of D L = 1225 Mpc (z = 0.252, 19]), and the mass 
within the scale radius of R « 800 kpc « 4.2arcmin is 
Mfgy 35 = 6.5 x 1O 14 M HJ. No other obvious mass 
concentrations are seen in the field of view. (Here, and 
in the following, a cosmology with Q, m = 0.26, £!a = 0.74, 
h = 0.71 is assumed). 

The spectral resolution of a grating spectrometer de- 
creases proportionally to the angular extension of the 
source. For Abell 1835 the emission is dominated by 
photons from within the core radius, r c — 39 ± 3 kpc « 
lOarcsec \lQ so this characteristic radius was used to 
determine the spectral resolution [34| : 
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where m is the dispersion order. Grating spectra of 
Abell 1835 were extracted from the observation with id 
511 using CI AO 3.3 j^]. The first order deflections were 
combined into a single spectrum with a spectral bin size 
of 0.006 keV and 0.013 keV (at £ 7 = 1 keV) for HEG and 
MEG, respectively. The line broadening due to the veloc- 
ity dispersion of dark matter in Abell 1835 is negligible 
(w/cw 10~ 5 - 10" 4 ). 

With the spectral fitting package Sherpa |2l|, a ther- 
mal plasma model (MEKAL, j23]) was fitted to the data 
with the temperature and abundance as free parameters. 
The model is an adequate representation of the data, giv- 
ing a reduced x 2 of 1.2 in both cases (683 and 569 degrees 
of freedom for HEG and MEG, respectively). The flux 
was determined using the slice method for a slice width 
given by eq. ^ 



FIG. 1: (Color online) Abell 520 observed in X-rays (0.3- 
10.0 keV, blue color) with Chandra with the gravitational 
potential from weak lensing overlaid (green contours). The 
dark matter blob in the red circle has very low X-ray emission 
from baryons. 



IV. THE DARK MATTER BLOB IN THE 
GALAXY CLUSTER ABELL 520 

For direct imaging data the field of view can be opti- 
mized by observing a dark matter dense region with low 
X-ray emission from baryons. A unique example is the 
merging cluster of galaxies Abell 520 containing a "blob" 
of high mass concentration with very low X-ray emission 
discovered recently using weak gravitational lensing |23| , 
see FIG. n]|3(l. 

An ACIS-S3 0.3-9.0 keV spectrum was extracted 
from the Chandra observation with id 4215 for 
a region centred at the dark blob (RA, DEC = 
04 h 54 min 04.036 s ,+02°52'37.30") and with a radius of 
r = 0.85arcmin = 190 kpc (the red circle in FIG. QJ. 
A model consisting of a power law and six Gaussians was 
fitted to the spectrum with a reduced x 2 °f 1-2 (for 80 
degrees of freedom). The flux was determined using the 
slice method for a slice width given by the resolution of 
ACIS-S3 El: 



a S3 = 0.005£ 7 + 0.05keV 



(2) 



The mass of the dark matter blob in A520 has been 
derived from weak gravitational lensing to be M^ l °^ = 
4.78±1.5x 10 13 ft. _1 M Q . This value is based on measuring 
the overdensity in the blob region with respect to the 
mean density in a surrounding annulus with inner and 
outer radius of 0.85 arcmin and 4 arcmin, respectively. 
Hence, the mass value can be regarded as a conservative 
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lower limit on the mass contained within the blob region. 
A detailed description of the data and methodology of 
the weak lensing analysis is given elsewhere [24^ . 

The luminosity distance to A520 is D L = 980 Mpc (z = 
0.203, H|). 



V. DECAY RATE 

Let us now consider a specific dark matter candidate, 
the sterile neutrinos. They can decay radiatively as v s — > 
where v a is an active neutrino. This is a two-body 
decay with a photon energy of E 7 = m s /2. Assuming 
only one kind of dark matter, the observed flux, F b Sl at 
a given photon energy yields an upper limit on the flux 
from two-body radiatively decaying dark matter: 
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FIG. 2: (Color online) The upper limit on the radiative 
two-body decay rate obtained from the dark matter blob of 
Abell 520 (dashed) and the grating spectra of Abell 1835 
(solid) shown together with the Milky Way halo constraint 
(ltH (dot-dashed). 



in Abell 520 it is M^ 20 



The Milky Way dark matter halo will always be in- 
cluded in the observation, and its mass contribution to 
the total mass in the field of view has to be taken into 
account [HI3. For the Milky Way we use a virial mass 
of M%% = 10 12 M and by integration of a NFW pro- 
file we find the halo mass to have a mean distance of 
35 kpc. The halo mass and mean distance varies less 
than a factor of two for a reasonable range of model pa- 
rameters 0| • For the grating observation of Abell 1835 
the Milky Way halo mass within the field of view is 
9 x 10 5 M Q , and for the dark matter blob 
4 x 10 3 M Q . 

FIG. |2] shows the upper limit on the decay rate of any 
dark matter candidate (given by eq. |HJ) obtained from 
the total amount of received flux for the Abell 520 dark 
matter blob and the Abell 1835 grating data. It is seen 
that the grating data of Abell 1835 provide constraints 
that are one to two orders of magnitude stronger than 
the constraints obtained from observations of the Milky 
Way halo alone [HIS- Also, the data from the blob 
of Abell 520 provide stronger constraints than the Milky 
Way halo. 

The constraints on the decay rate from the grating 
data of Abell 1835 can be approximated by an analytical 
expression for HEG and MEG independently. As seen 
in FIG. a second order polynomial describes the HEG 
constraints quite well, and a sixth order polynomial fits 
the MEG data: 
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FIG. 3: (Color online) The Abell 1835 grating spectrum up- 
per limit on the decay rate and its approximate analytical 
expression given by of eq. 0] 



VI. CONSTRAINING MASS AND MIXING 
ANGLE 

By regarding all decay branches possible through os- 
cillations, the mean lifetime of a sterile Dirac neutrino of 
m,, has been determined to be I26ll27l: 



mass, 



f(m s ) ■ 10 



20 



Ttot (m s /kev) sin 2 (26») 



(5) 



where Ttot is the total decay rate. f{m s ) takes into ac- 
count the open decay channels so that for m s < 1 MeV, 
where only the neutrino channel is open, f(m s ) = 0.86. 
For Majorana neutrinos, which we will be considering 
below, f(m s ) is half the value for Dirac neutrinos. The 
branching ratio of the radiative decay is r 7 /r to4 = 
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FIG. 4: (Color online) The observational constraints from 
the grating spectrum of Abell f 835 (red) and the dark matter 
blob of A520 (orange) together with the Tremaine-Gunn limit 
(black, |2Sjp and earlier X-ray constraints (yellow, [l3l . ll5l . H6l 

111 Hi). 

27a/8?r w 1/128 [H. This can be combined with eq. 
to give: 

sin 2 (20) < TxlO^sec- 1 ( f d %. ) (^)"(6) 
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The observational constraints in the sin 2 (28) — m s 
parameter space are shown in FIG. 0] for the grat- 
ing spectrum of Abell 1835 (red) and the dark matter 
blob of A520 (orange) together with the Tremaine-Gunn 
limit (blac k. 1281) and earlier X-ray constraints (yellow, 

[3 HE Ella), it is seen that even though no op- 
timization of the field of view can be performed for the 
grating data, the improvement of instrumental spectral 
resolution leads to superior constraints. The grating con- 
straints are very robust as they have been derived from 
the total amount of received X-ray flux without subtrac- 



tion of any background contributions. The earlier con- 
straints are derived by less robust and more model de- 
pendent methods. 



VII. CONCLUSIONS 

A very general constraint on the decay rate for all dark 
matter particle candidates with a two-body radiative de- 
cay in the X-ray range has been derived. We have an- 
alyzed grating spectra of the galaxy cluster Abell 1835, 
and a spectrum obtained through direct imaging of the 
almost pure dark matter blob in the galaxy cluster 
Abell 520. These provide the strongest constraints on 
radiatively decaying dark matter derived in a very conser- 
vative way from cosmic X-ray data. The mass and mixing 
angle can be constrained in the specific case of sterile neu- 
trinos, leaving a low mass (m s < 10 keV) and low mixing 
angle (sin 2 (28) < 10 -6 ) window open. However, com- 
bined with other methods [3(], OH) , the Dodelson-Widrow 
scenario for sterile neutrino production could soon be 
ruled out. Other mechanisms for production of sterile 
neutrinos, for example, via inflaton decays [3^, would 
satisfy all the constraints. 

The obtained constraints can be improved significantly 
by improving the signal to noise ratio (optimization of 
field of view) and by improving the instrumental spectral 
resolution. This could be realized through X-ray grating 
observations of dark matter dense regions or the Milky 
Way halo (blank sky fields). 
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